function c = fn_c(w,tau,n,g)

c = (1-tau).*w.*n + g;

end